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ABSTRACT 


A  procedure  is  proposed  for  analysing  electromagnetic  interference  (EMI)  with 
Middleton  models.  An  introduction  to  the  models  is  given,  and  descriptions 
of  automated  procedures  for  estimation  and  optimisation  of  model  parameters 
are  proposed.  Various  techniques  are  developed  to  solve  numerical  problems 
in  the  calculation  of  the  model  parameters.  Many  of  the  practical  problems 
associated  in  the  analysis  are  outlined. 
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A  Practical  Utilisation  of  Middleton  EMI  Models: 
Automated  Modelling,  Parameter  Estimation  and 

Optimisation 

EXECUTIVE  SUMMARY 

This  report  presents  the  results  of  an  analysis  of  high-frequency  electromagnetic  inter¬ 
ference  (HF  EMI  or  ‘radio  noise’).  The  aim  of  the  analysis  was  to  develop  software  to 
model  HF  EMI  and  extract  useful  physical  information  from  those  models  as  a  basis  for 
simulation.  Sections  concerning  data  acquisition  and  processing  are  present. 

The  software  successfully  modelled  the  HF  EMI  using  Middleton’s  models.  However, 
the  extraction  of  useful  physical  information  was  limited,  as  the  models  are  quite  com¬ 
plex  and  information  available  in  the  literature  is  incomplete.  The  modelling  procedures 
nevertheless  provide  a  useful  insight  into  HF  EMI  and  could  form  the  basis  for  modem  per¬ 
formance  simulation  tests,  using  information  about  interferer  positions,  densities,  trans¬ 
mitting  powers,  etc.  That  is,  interference  would  be  simulated  entirely  by  approximate 
physical  situations,  compared  with  traditional  methods  using  physically  meaningless  pa¬ 
rameters. 
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1  Introduction 

Middleton’s  EMI  models  aim  to  provide  realistic,  quantitative  descriptions  of  man¬ 
made  and  natural  EMI  environments,  accurate  even  with  highly  non- Gaussian  random 
processes.  All  model  parameters  can  be  extracted  from  measurement,  are  physically 
meaningful,  and  the  model’s  structures  are  based  on  underlying  physics  of  radio  wave 
propagation  and  reception.  Hence  knowledge  of  the  model’s  parameters  provides  insight 
as  to  the  physical  processes  taking  place.  What  this  in  effect  means  is  that  the  EMI  process 
resulting  from  the  sources  incident  on  the  receiver  may  be  simulated  with  the  knowledge 
of  several  key  parameters. 

The  authors  aim  firstly  to  impart  a  basic  knowledge  of  the  Middleton  EMI  models 
(omitting  the  mathematical  derivations)  and  to  explain  in  exactly  what  way  its  parameters 
are  physically  meaningful.  Methods  are  then  proposed  to  enable  extraction  of  estimates 
of  model  parameters  from  data  samples,  and  to  optimise  these  parameters.  The  question 
“How  well  does  the  model  fit  the  empirical  data”  is  answered.  Many  of  the  numerical 
problems  and  solutions  associated  with  such  a  process  are  given,  however,  details  of  the 
actual  computer  algorithm  structures  are  omitted.  The  paper  also  explains  in  simple  terms 
how  the  algorithms  decide  on  which  model  (Class  A  or  B)  to  use,  and  how  this  affects  the 
estimation  and  optimisation  routines. 

There  are  several  other  non-Gaussian  EMI  models  that  may  be  of  use.  The  Kabanov 
[13,  pp,  122-129]  and  K-distributed  [14]  models  have  shown  promise:  what  may  be  con¬ 
sidered  a  flaw  in  these  two  models  is  their  lack  of  insight  into  the  physical  mechanisms 
behind  the  EMI.  It  is  for  this  reason  that  Middleton  models  are  favoured  for  analysis. 


2  The  Middleton  EMI  Models 

2.1  General  Description 

There  are  three  classes  of  Middleton  models:  A,  B,  and  C.  Class  A  models  are  used 
when  the  noise  bandwidth  is  narrower  than  the  receiver  bandwidth,  producing  ignorable 
transients  in  the  receiver’s  front-end  RF  and  IF  (radio  and  intermediate  frequency  respec¬ 
tively)  stages.  Class  B  models  are  used  when  the  incoming  noise  is  spectrally  wider  than 
the  receiver  bandwidth,  so  that  transient  effects  in  build  up  and  delay  occur  to  produce 
‘ringing.’  Class  C  EMI  may  occur  when  both  Class  A  and  B  interferers  are  present.  How¬ 
ever,  as  information  available  on  Class  C  models  is  severely  limited,  comments  concerning 
them  shall  not  occur  again  in  this  discussion. 

The  Middleton  models  are  supposedly  effective  for  naturally  occurring  or  man-made 
noise  processes  (occurring  separately  or  simultaneously)  for  signals  containing  intelligent 
or  random  components.  They  fail  when  the  data  samples  are  either  (i)  nonstationary1,  or 
(ii)  inhomogeneous. 

1In  the  context  of  this  paper,  stationarity  is  assumed  to  be  present  if  the  underlying  probabilistic 
mechanisms  are  time-invariant  during  the  data  acquisition  period.  For  example,  there  are  no  changes  in 
average  emitter  numbers,  propagation  laws,  emitter  locations,  etc.  A  change  in  this  mechanism  during 
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Worthy  of  note  is  the  fact  that  the  models  assume  incident  signal  phases  to  be  uni¬ 
formly  distributed  on  [0,  2vr].  With  multipath  effects,  this  may  not  be  the  case.  Also,  the 
models  assume  two  independent  components:  (i),  Poisson-distributed  (in  space  and  time) 
‘impulsive’  EMI,  containing  only  a  few  (ten  or  less)  discrete  sources  of  relatively  high  level, 
and  (ii),  zero-mean  Gaussian  background  noise  from  either  receiver  (self-generated)  noise 
or  external  noise.  There  are  five  more  assumptions  regarding  the  physical  processes  of 
source  emissions;  however  these  may  prove  less  crucial  (refer  [3,  p.8.])  They  are;  (a)  a 
far-field  (Fraunhofer)  condition  exists  [8];  (b)  there  is  little  relative  motion  (viz.  Doppler 
shift)  between  sources  and  receiver;  (c)  the  typical  source  has  a  beam  pattern  that  is  not 
necessarily  omni-directional;  (d)  the  receiving  antenna  generally  has  a  directional  beam 
pattern;  and  (e)  the  source  may  have  a  time- variable  probabalistic  mechanism  (e.g.,  change 
in  level,  frequency). 

Because  of  the  form  of  the  EMI,  it  becomes  useful  to  graphically  represent  the  data  on 
what  is  called  a  Rayleigh  Probability  Presentation  (RPP),  which  is  basically  a  log-log  by 
log  scale  (specifically,  it  is  a  logio[-ln[...]]  vs  logio[...]  scale)  plot.  The  RPP  shows  many 
characteristics  of  the  signal’s  envelope,  including  the  percentage  of  time  it  is  Rayleigh- 
behaved  (as  is  thermal  noise).  Thermal  (Gaussian,  or  Rayleigh-enveloped)  noise  appears 
as  a  straight  line  with  gradient  on  a  RPP.  The  vertical  axis  displays  the  quotient  of 
the  envelope  voltage  E  and  the  root  mean  square  of  the  envelope  voltage  ETma,  while  the 
horizontal  axis  displays  the  percentage  of  time  a  certain  envelope  threshold  E0  is  exceeded, 
and  is  hence  referred  to  as  ‘exceedance  probability’.  The  plot  is  therefore  known  as  an 
Amplitude  Probability  Distribution  (APD). 


2.2  The  Middleton  Class  B  EMI  model 

Figure  1  shows  a  typical  APD  for  Class  B  EMI.  Shown  clearly  are  two  regions,  BI 
and  BII,  separated  by  an  inflection  point  IP.  e  and  6o  are  the  normalised  envelope  voltage 
and  the  corresponding  threshold  respectively.  The  composite  model  may  be  considered  an 
approximation  for  low  and  high  values  of  the  envelope.  Region  I,  which  is  an  approximation 
for  low  envelope  values  is  described  by  the  equation: 


Pj5/[e  >  €0]  =  1  —  eo2  X  ^ 


(-1)"  x  A2  X  r(l  +  ^)  x  M(1  +  ^;2;-6q2) 


n=0 
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eo  <  es  (1) 
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data- acquisition  would  be  reflected  by  a  changing  APD  and  hence  changing  model  parameters. 
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Eo/Erms 


Figure  1:  Typical  models  (Class  B  model  parameters:  a=1.8,  Qb=0-47,  $b=5,  Nj=1.63, 
Aa=l,  Ab  =0. 0001,  Class  A  model  parameters:  Aa=0.001,  Qa=1,  ^a—1  ) 


eo  -  ( 


go  x  Nt 
2xG  ’ 


and  eo 


-  ...  =  =  normalised  envelope  threshold 

y/2  X  «B  X  (1  +  §B) 


a  =  spatial  density  propagation  parameter 
0.B  =  intensity  of  the  impulsive  component 

$B  =  ratio  of  the  intensities  of  the  Gaussian  to  the  non-Gaussian  components 
Eq  =  un-normalised  envelope  threshold 
Nj  =  scaling  parameter 
T  =  Gamma  function 

M  =  Rummer’s  confluent  hypergeometric  function 
Aa  =  scaling  parameter 

Region  I  has  a  finite  non-zero  gradient  on  the  APD  at  low  probabilities,  and  as  such 
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does  not  satisfy  the  condition  for  finite  energy  (and  hence  the  existence  of  all  finite  moments 
of  the  envelope.)  The  required  ‘bend-over’,  provided  by  Region  II,  has  the  form: 


2.3  The  Middleton  Class  A  EMI  model 

Class  A  models  are  relatively  more  straight-forward  to  derive  than  Class  B  models2. 
Visually,  it  is  simple  to  distinguish  between  the  two  on  an  APD,  Class  A  having  a  more 
pronounced  ‘elbow’  region  than  Class  B.  An  example  is  shown  in  Figure  1.  It  is  not  a 
composite  model:  the  one  equation  describing  it  being; 


Pa[c  >  eo] 


(3) 


where: 


(**  +  £) 

(1  +  *a)  ’ 


where  and  Aa  are  analogous  to  their  respective  Class  B  parameters. 

Class  B  models  will  usually  work  well  with  atmospheric  and  ignition  noise,  etc.,  while 
(in  general)  Class  A  will  work  with  ‘intelligent’  (that  is,  man-made  narrowband  EMI 
meant  to  convey  information — eg,  modems.)  noise  processes.  This  depends  on  the  receiver 
bandwidth,  which  is  usually  assumed  wider  than  narrow-band  intelligent  EMI. 

2  A  point  worth  noting  is  that  Class  A  EMI  is  a  special  case  of  (Region  II)  Class  B  EMI,  as  a  -v  oo. 
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2.4  Physical  Significance  of  Model  Parameters 

The  physical  meaning  of  these  parameters  may  now  be  given,  beginning  with  Region 
I  in  the  Class  B  model.  $a,B  is  the  ratio  of  the  intensities  of  the  Gaussian  to  the  non- 
Gaussian  components3;  on  the  APD  it  is  the  Eo/Erms  value  where  the  data  begins  to 
depart  from  Rayleigh  (straight-line)  behaviour.  The  actual  values  of  the  above  intensities 
given  above  may  be  found  quite  easily,  by  methods  explained  below.  Another  parameter, 
fi A,b ,  that  is  not  shown  above  but  is  used  in  the  process,  is  the  intensity  of  the  impulsive 
component.  The  mean-square  envelope  of  the  resultant  signal  can  be  calculated  by  using 
the  above  two  (shape)  parameters: 


E2  =  20^)b(1  +  $a,b) 


Aa,b  is  the  impulsive  index:  the  product  of  the  average  number  of  source  emission 
events  detected  by  the  receiver  and  the  mean  duration  of  a  typical  intefering  source  emis¬ 
sion.  Very  large  values  indicate  a  high  density  of  overlapping  waveforms  at  a  given  instant, 
giving  an  (asymptotically)  Gaussian  form  to  the  EMI,  as  is  shown  in  Figure  2(a).  Low 
values  indicate  little  overlapping,  and  hence  the  contributions  of  only  a  few  sources  are 
significant  at  a  given  instant,  as  is  shown  in  Figure  2(b).  This  case  would  imply  a  highly 
impulsive  resultant  waveform.  Graphically,  Aa,b  is  the  exceedance  probability  value  where 
the  inflection  occurs,  and  as  such  this  property  of  the  signal  may  be  easily  seen  at  a  glance. 

a  (0  <  a  <  2)  is  the  spatial  density  propagation  parameter,  since  its  value  depends 
on  the  interacting  spatial  effects  of  source  density  and  source  propagation.  That  is,  the 
effective  source  density  and  the  propagation  law  (signal  attenuation  with  distance)  are 
both  functions  of  a.  Graphically,  it  is  a  complex  measure  of  the  gradient  of  the  empirical 
APD  in  a  small  interval  below  the  inflection  point.  Point  A  is  used  as  a  source  in  this 
example,  and  is  d  units  from  the  receiver,  a  allows  the  following  to  be  said  about  the 
receiver  power  Pr: 


Pr  oc 


1 

d*-“ 


where  i  —2  or  3,  depending  on  the  source-receiver  path  medium. 

Hence  a  gives  the  power-law  exponent  used  for  the  particular  propagation  medium. 
Also,  the  value  of  a  depends  upon  the  source  density.  However,  it  is  not  clear  how  the 
source  domain  is  determined. 

Ni  and  Aa  are  scaling  parameters,  required  in  joining  the  two  regions  to  form  a  smooth 
composite  pdf.  As  such,  they  convey  no  physical  meaning.  Decreasing  the  Nj  value  will 
raise  the  region  I  Eo/Erms  value  without  affecting  the  exceedance  probability  value. 
Increasing  the  Aa  value  will  also  increase  the  Region  I  Eo/Erms  value,  but  will  affect  the 
exceedance  probability  value. 

3A  parameter  with  two  subscripts  indicates  relevance  to  both  Class  A  and  B  models. 
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Figure  2:  Source-receiver  geometry,  •-sources,  receiver . 


For  some  idea  of  how  the  parameters  affect  the  shape  of  the  model  on  the  RPP,  refer 
to  [5,  Figures  2. 1-2.4,  3.1,  and  3.2]  and  especially  [4,  pages  91-99]4. 

3  Parameter  Estimation 

The  parameters  are  estimated  according  to  their  graphical  significance  by  Middle¬ 
ton’s  method  [6,  especially  pages  197-201].  For  example,  $AtB  is  simply  found  by  noting 
the  Eo/Erms  value  where  the  data  departs  from  Rayleigh  behaviour.  This  is  relatively 
straight-forward,  as  are  the  methods  of  finding  the  other  key  parameters.  Note  that  al- 

4  A  more  comprehensive  study  of  the  effects  of  parameter  variation  is  given  in  the  authors’  notes  on  the 
subject. 
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though  visually  it  may  be  easy  to  estimate  the  parameters,  the  algorithm’s  automated 
processes  need  to  eliminate  all  subjectivity  involved  in  such  a  method,  which  can  become 
a  problem  with  some  of  the  more  unusual  data  samples. 

Parameter  estimation  for  Class  A  is  quite  easy,  as  there  are  only  three  parameters.  §  a 
is  estimated  by  the  method  mentioned  above;  A  a  by  finding  the  inflection  point,  is 
then  scaled  until  the  model  and  data  appear  more  closely  matched  on  the  APD. 

Class  B  parameter  estimation  can  become  difficult  at  times  because  of  its  complexity 
(six  parameters).  §b  and  Ab  are  found  by  a  similar  method  as  above,  however  the 
remaining  parameters  introduce  difficulty,  a  and  f Ib  are  given  an  arbitrary  (but  sensible; 
in  the  middle  of  their  allowable  ranges)  value  as  a  starting  point.  These  four  parameters 
describe  Region  II,  so  they  are  passed  into  an  optimisation  routine  specifically  for  this 
region.  Once  these  four  parameters  are  found  it  becomes  necessary  to  estimate  the  values 
of  Ni  and  Aa  that  will  result  in  a  smooth  composite  function  around  the  inflection  point. 
This  is  done  by  forcing  two  requirements  upon  the  model: 

(!)  in  the  Rayleigh  region. 

(2)  Pbi  —  Pbii  at  the  inflection  point. 

Requirement  (1)  states  that  in  the  straight-line  Rayleigh  section  of  the  RPP,  the  two 
regions  must  have  the  same  gradient,  and  requirement  (2)  stipulates  that  the  two  regions 
must  have  the  same  value  at  the  composite  function’s  inflection  point.  Applying  the  above 
arguments  to  Equations  1  and  2,  it  can  be  shown  that: 


with  Aa  as  in  Equation  1,  and  as  as  in  Equation  2. 

Aa  is  found  by  solving  for  the  root  of  Equation  4.  This  form  cannot  be  solved  al¬ 
gebraically,  so  an  iterative  procedure  was  written  to  find  the  numerical  value.  Once  the 
value  of  Aa  is  found,  Ni  may  be  found  easily  by  use  of  the  following  equation  derived  from 
requirement  (1); 


7 


DSTO-TR-0234 


Ar  e-^  X  fa  ~  1 

Kl  = - A - 


(5) 


This  completes  the  parameter  estimation  and  allows  the  formulation  of  a  smooth  com¬ 
posite  function.  However,  the  above  method  for  estimating  Aa  and  JVj  is  not  accurate 
enough;  another  optimisation  routine  is  needed  to  optimise  these  two  parameters,  using 
the  remaining  four  parameters  to  implement  the  Region  I  part  of  the  model. 


The  above  procedure  is  one  suggested  method  for  parameter  estimation.  Others  exist, 
some  estimating  on  a  sample- by- sample  basis  [16]. 


4  Model  Implementation 

The  model  is  implemented  by  various  routines  using  the  “Mathematica”  software  pack¬ 
age.  The  routines  perform  the  following:  (i)  estimation  of  the  parameters,  (ii)  modelling  of 
the  data,  (iii)  a  (limited)  optimisation  of  the  parameters,  and  (iv)  characterisation  of  the 
data  sets.  The  most  crucial  step  is  selecting  the  inflection  point.  Knowledge  of  this  allows 
a  graphical  categorisation  of  the  APD  as  Class  A  or  B;  a  large  gradient  at  the  inflection 
point  implying  Class  A,  while  a  smaller  gradient  implying  Class  B.  In  reality,  both  models 
are  fitted  to  the  data  because  of  a  number  of  borderline  cases:  it  becomes  difficult  to  label 
some  either  way.  Class  B  models  are  orders  of  magnitude  slower  to  evaluate  than  Class  A 
because  of  the  need  to  evaluate  Gamma  and  Hypergeometric  functions. 

The  routines  test  the  validity  of  the  model  with  respect  to  the  data  sets.  This  involves 
a  95%  significance  level  Kolmogorov- Smirnoff  (K-S)  goodness-of-fit  test  [2],  which  is  rec¬ 
ommended  by  Middleton  [7,  page  214.]  Other  procedures  test  whether  either  model  will 
be  valid  under  all  circumstances,  by  (amongst  others)  testing  the  angle  of  the  inflection 
point5. 


5  Parameter  Optimisation 

Optimisation  (for  our  purposes)  is  the  minimisation  of  the  sum  of  the  squares  of  the 
errors  between  the  dataset  APD  and  the  model  APD.  Hence  optimisation  of  the  model 
parameters  shapes  the  model  APD  in  order  to  minimise  this  statistic. 

Fortran  77  programs  are  used  to  perform  the  optimisation  using  quasi-Newton  and 
finite- difference  gradient  methods  [10].  These  programs  operate  at  much  faster  rates  than 
those  written  in  Mathematica,  with  the  disadvantage  that  several  numerical  evaluation 

B  A  good  measure  of  whether  a  Class  B  model  will  fit  a  data  set  is  the  maximum  (transformed)  difference 
of  the  exceedance  probabilities  between  the  data  and  an  optimised  Class  A  model.  If  it  is  above  a  certain 
threshold,  a  Class  B  model  will  not  be  possible,  and  hence  one  or  more  assumptions  are  erroneous. 
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problems  present  themselves6. 

The  optimisation  of  Class  A  model  parameters  is  quite  straight  forward  and  does  not 
require  special  optimisation  software.  Region  II  of  the  Class  B  model,  although  having 
an  extra  parameter  cr,  behaves  in  exactly  the  same  way  as  a  Class  A  model.  The  extra 
parameter  still  allows  simple  and  predictable  adjustment  of  the  APD  in  this  Region. 

The  principle  difficulty  lies  in  optimising  Region  I  of  the  Class  B  model.  The  shape 
and  scale  effects  of  the  four  parameters  here  are  not  well  understood  or  predictable,  and 
specially-written  optimisation  software  must  be  used.  It  was  found  that  the  initial  es¬ 
timates  of  the  parameters  need  to  be  carefully  selected  for  the  optimisation  to  converge 
correctly.  Haphazard  selection  of  the  estimates  may  lead  to  an  optimisation  routine  con¬ 
verging  to  a  local  (rather  than  global)  minimum. 

Upon  optimising  the  model,  the  implementation  of  the  figure  of  merit  for  the  goodness 
of  fit  should  be  thought  of  carefully.  The  K-S  test  is  insensitive  to  rare-event  discrepancies 
(difference  between  data  and  model)  the  less  ‘rare’  they  become.  The  rare  events  are 
important  in  the  process,  because  of  the  need  to  include  impulsive  behaviour  in  the  models, 
and  hence  a  K-S  test  on  the  whole  exceedance  probability  range  is  inappropriate.  What 
is  needed  is  a  statistical  test  that  weights  increasingly  smaller  probabilities  more  heavily. 
The  Anderson-Darling  statistic  has  this  property  and  may  be  effective  [15,  Sec.  VIII]. 


6  Data  Acquisition 


Data  acquisition  is  a  vital  step  before  the  modelling  process,  and  issues  such  as  sampling 
rates  and  stationarity  must  be  addressed  and  resolved  before  modelling  takes  place. 

Ideally,  the  sampling  rate  should  be  above  the  Nyquist  rate  (twice  the  highest  frequency 
component)  to  avoid  aliasing  [11,  1.4.1  and  1.4.2].  Remembering  that  the  Class  B  model 
is  a  model  of  band-limited  EMI  anyway,  it  becomes  apparent  that  this  does  not  pose 
a  problem.  However,  the  sampling  (and  hence  filtering)  should  be  performed  such  that 
the  high-frequency  impulsive  nature  of  the  EMI  is  not  drastically  altered.  As  is  shown 
elsewhere  [12,  pp.  44,45],  excessively  narrow  filter  bandwidths  cause  attenuation  of  rare 
events  (say  those  around  0.0001%  exceedance  probability  on  the  APD)  by  as  much  as  20 

6Primarily,  the  problem  is  the  implementation  of  Kummer’s  Confluent  Hypergeometric  function  from 
first  principles.  Various  solutions  are  given  in  text  [9],  however  the  authors  have  derived  a  recursive 
equation,  shown  below,  that  eliminates  several  numerical  evaluation  problems: 


OO 

M(a,  b,  z)  =  ^  Mn , 

n=0 

where: 


Mn  = 


1 


az 


(q-fn  —  l)z 
(6+n  — l)n 


M„_ 


l 


if  n  =  0 
if  n  —  1 
if  n  >  1 
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dB  by  changing  the  bandwidth  by  only  several  kHz.  A  change  of  this  magnitude  may  well 
render  a  Middleton  model  ineffective.  The  attenuation  will  take  the  form  of  high,  narrow 
power  peaks  being  amplitude-clipped  and  pulse- width  extended.  Hence  persistance  of ‘flat’ 
APD’s  of  apparently  Class  B  form  may  be  caused  by  low  sampling  rates. 

It  was  suspected  that  autocorrelation  caused  by  oversampling  would  affect  the  prob¬ 
ability  values  on  the  APD.  To  study  the  effects  of  autocorrelation,  every  nth  value  was 
taken,  where  n  ranged  from  2  to  10.  It  was  found  that  the  effects  were  negligible:  taking 
the  nth  sample  did  not  appreciably  alter  the  APD  below  around  n  =  10,  and  hence  it  was 
assumed  autocorrelation  did  not  pose  a  problem  in  the  analysis. 

It  is  important  to  have  stationary  data  for  the  modelling  process,  as  it  is  one  of  the 
assumptions  Middleton  makes.  Nonstationarity  becomes  more  of  a  problem  the  longer 
the  data  acquisition  period  is,  however,  this  (it  would  seem  excessive)  time  is  needed  to 
establish  the  ‘rare  event’  portion  of  the  sample  set.  This  is  because  with  large  (but  finite) 
sample  sets,  even  in  the  order  of  hundreds  of  millions  of  samples,  the  ‘rare  events’  portion 
of  the  model  will  not  necessarily  represent  the  respective  portion  of  the  EMI  closely.  There 
may  in  the  order  of  ten  or  less  samples  representing  the  data  in  these  portions7.  It  would 
seem  that  a  compromise  must  be  reached  in  this  process:  dividing  the  data  into  smaller 
intervals  giving  different  APD’s  implies  nonstationarity  of  the  whole  interval,  but  one  can 
only  divide  so  far,  as  APD’s  can  only  accurately  represent  the  data  when  they  are  above 
a  certain  sample  size. 

Specific  examples  of  modelling  may  be  seen  in  another  report  [1].  This  report  shows 
in  more  detail  the  practical  problems  associated  with  modelling  and  simulation. 

7  Conclusions  and  Recommendations 

The  procedures  outlined  in  this  paper  provide  the  basis  for  an  elementary  analysis 
of  the  APD’s  of  EMI.  Middleton’s  papers  alone  may  not  do  so.  The  bulk  of  his  papers 
give  intricate  details  as  to  the  mathematics  of  devising  the  various  models:  they  do  not 
give  a  good  foundation  for  a  practical  utilisation  of  the  models,  nor  do  they  approach 
the  problems  associated  with  such  a  task.  As  Middleton’s  papers  are  both  repetitious  and 
confusing,  the  authors  have  attempted  to  explain  the  models  in  simpler  terms,  disregarding 
most  theory  and  focusing  on  explanation  and  implementation. 

The  Middleton  model  approach  offers  substantial  advantages  over  other  analysis  and 
characterisation  methods;  the  primary  one  being  the  physical  nature  of  the  models.  With 
a  good  understanding  of  the  theory  and  implementation,  Middleton  modelling  may  prove 
to  be  a  valuable  tool. 


7For  example,  sampling  at  10  kHz  for  10  minutes  will  produce  an  APD  with  the  rare  events  portion 
(say  0.0001%)  being  represented  by  (on  average)  only  6  samples. 
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